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ABSTRACT 

We analyze the mean rest-frame ultraviolet (UV) spectrum of Type la Supernovae (SNe) and its 
dispersion using high signal-to- noise Keck-I/LRIS-B spectroscopy for a sample of 36 events at in- 
termediate redshift (z=0.5) discovered by the Canada-France-Hawaii Telescope Supernova Legacy 
Survey (SNLS). We introduce a new method for removing host galaxy contamination in our spectra, 
exploiting the comprehensive photometric coverage of the SNLS SNe and their host galaxies, thereby 
providing the first quantitative view of the UV spectral properties of a large sample of distant SNe la. 
Although the mean SN la spectrum has not evolved significantly over the past 40% of cosmic history, 
precise evolutionary constraints are limited by the absence of a comparable sample of high quality 
local spectra. The mean UV spectrum of our z ~0.5 SNe la and its dispersion is tabulated for use in 
future applications. Within the high-redshift sample, we discover significant UV spectral variations 
and exclude dust extinction as the primary cause by examining trends with the optical SN color. 
Although progenitor metallicity may drive some of these trends, the variations we see are much larger 
than predicted in recent models and do not follow expected patterns. An interesting new result is a 
variation seen in the wavelength of selected UV features with phase. We also demonstrate systematic 
differences in the SN la spectral features with SN lightcurve width in both the UV and the optical. 
We show that these intrinsic variations could represent a statistical limitation in the future use of 
high-redshift SNe la for precision cosmology. We conclude that further detailed studies are needed, 
both locally and at moderate redshift where the rest-frame UV can be studied precisely, in order that 
future missions can confidently be planned to fully exploit SNe la as cosmological probes. 
Subject headings: surveys - supernovae: general ~ cosmological parameters 



1. INTRODUCTION 

Supernovae of Type la (SNe la) are now well- 
established as cosmological distance indicators. In ad- 
dition to the original surveys by the Supernova Cos- 
mology Project (SCP: lPerlmutter et al.llT9 97. 199 g) and 
the High-Z Supernova Search Team (|Schinidt et al.[l998l : 
[Riess ct al. 1998), a new gener ation of SN la sur- 
veys is underw a y bot h locallv (|Aldering et al.l 120021 : 
iLi &: Filippenkol 120051: iHamuy et al I I2006D an d at 
higher redshifts ([Astier et al.l 120061: iRiess et al] l2007t 
IWood-Vasev et all l2007| ). Despite the availability of 
independent probes of the presence and properties of 
dark en ergy from studies of the cosmic microwave back- 



xound (ISpergel et al.l 2007f) and galaxy redshift surveys 
Efstathiou et al.ll2002l : ICole et aLlbooa lEisenstein et all 
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l2005f ). the luminosity distance-redshift relation for 
SNe la provides the only direct evidence for a cosmic 
acceleration. 

The detection and characterization of dark energy, via 
measurements of the average cosmic equation of state 
parameter <w>, requires the precision measurement of 
SNe la to redshif ts z :^0. 5-1, samp ling the epoch of cos- 
mic acceleration (jAstier et al.ll2006[) . However, more pre- 
cise constraints on the nature of dark energy, for example 
evidence for any variation in w with reds hift, requires ex- 
tendi ng these studies to redshift z>l (jRiess et al.ll2004 
l2007f ) where the early effects of deceleration may be de- 
tectable. As projects are developed whic h plan to probe 
SNe la beyond z=l for this purpose (e.g lAlderin i [2005l : 
iBenford fc Lauei1l2006f ). it becomes important to under- 
stand the possible limitations of using SNe la as distance 
probes. Key issues relating to the diversity of SNe la and 
their possible evolution with redshift as a population, to- 
gether with the limiting effects of dust and/or color cor- 
rections to their photometric properties, are particularly 
crucial to understand. 
Sev eral local studies dHamuv et al.l 119951: iRiess et al 



T^, 'Hamuv et al.' '2000': ' HoweU 120011: iMannucci et ~ 



2005; Gallagher et al. 2005) have already indicated cor- 
relations between SN la properti es and ho st galaxy mor- 
phologies. More recently, Sullivan et al.l ([2006b) have 
shown that the properties of distant SNe la appear to 
be a direct function of their local stellar population, 
with the distribution of light curve widths and hence 
peak luminosities correlating with the host galaxy spe- 
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cific star-formation rate. This work also determined that 
the rate of SNe la per unit stellar mass of their host 
galaxies is larger in actively star- forming galaxies, sug- 
gesting many must be produced quite rapidly in recently 
formed stellar populations, perhaps suggestive of more 
than one progenitor mechanism. The authors conclude 
that SNe la may well be a bim odal or a more complex 
populatio n of events (see also IScannapieco fc BildstenI 
[2005; Ma nnucci et al.|[2"006D . Such diversity in the prop- 
erties of SNe la could have far-reaching implications, par- 
ticularly if the mix of mechanisms or delay times within 
the broad population gradually changes with lo ok-back 
time (e.g. iMannucci et al.l l200d iSullivan et all l2006bl : 



IHowell et al.ll2007D . 

These recent developments, which illustrate how im- 
proved precision reveals new physical correlations in the 
SN la population, raise the broader question of whether 
future SN la experiments might be limited in precision by 
variations of a systematic nature within the population, 
for example with redshift, which cannot be removed via 
empirical corr elations. D etailed local surveys such as the 
LOSS/KAIT (iLi fc Filip pcnko 2005) and CfA surveys 
(jRiess et"al1ll999t iJha et al..,2006 ') have presented valu- 
able data on the homogeneity and trends in the SNe la 
population. Further promising work is being undertaken 
via the Supernova Factory ( Al dcring _ct_aL ,_2Q02|) a nd the 
Carnegie Supernova Project (jHamuv et al.l l2006l ) . Im- 
portant though these continued programs will be, they 
are insufhcient to address all possible concerns about the 
use of SNe la as precision tools in cosmology. Compara- 
ble studies at intermediate redshift^ will be particularly 
important in order to address questions relating to possi- 
ble evolutionary effects and environmental dependencies. 
In addition, it is not always practical at low redshift to 
cover the full wavelength range necessary to test for sys- 
tematic trends. 

In this paper we analyze high signal-to-noise ratio 
rest-frame ultraviolet (UV) spectra of a large sam- 
ple of intermediate redshift SNe la drawn from the 
Canada- France-Hawaii Telescope Supernova Legacy Sur- 
vey (SNLS; Astier et al. 2006), a rolling search which 
is particularly effective for locating and studying events 
prior to their maximum light. Our aim is obtain a 
substantially higher signal-to-noise in the spectra than 
that typically obtained during spectroscopic programs 
to type SNe and measure redshifts. We target the 
UV wavelength region because in this wavelength region 
the SN spectrum is thought to provide t he most sen- 
sitive probe of progeni tor metallicity (e.g. iHoflich et al.l 
119981 : iLentz et al.ll200Cil ) , a variable which may shed light 
on the possibility of progenitor evolution. The time- 
dependent UV spectrum is also needed in estimating 
"cross-band" fc-corrections, particularly at redshifts z > 
1 w here optical bandpasses probe the rest-frame near- 
UV (jRiess et al.ll2004 l2007f ). Little is known about the 
properties and homogeneity of the UV spectra of SNe la, 
largely because of the absence of suitable instruments for 
studying this wavelength range in local events. Although 
some local SN la UV spectra are available from Interna- 
tional Ultraviolet Explorer (lUE ) or Hubble Space Tele- 
scope (HST) satellite data (e.g. iLeibundgut et al.lfr991l : 
iKirshner et"aI1ll993t ICappellaro et al.lll995( l. the bulk of 

^ Defined here to represent the range 0.2<2<0.7. 



the progress now possible in this area can be provided 
from optical studies of intermediate-redshift events with 
large ground-based telescopes. 

The goals of this paper are thus to address the ques- 
tion of both the diversity and possible physical evolution 
in the intermediate redshift SN la family. We compare 
the rest-frame UV behavior of local SNe la with that de- 
rived for intermediate-redshift (z ~0.5) events where the 
rest-frame UV enters the region of high efficiency of the 
Keck LRIS-B spectrograph. We also study the degree to 
which the UV spectra at intermediate redshift represent 
a homogeneous population, independent of other vari- 
ables such as the physical environment and light curve 
stretch. 

A plan of the paper follows. In § [2] we introduce the 
salient features of the SNLS and our method for select- 
ing SNe la for detailed study. In §[3] we discuss the Keck 
spectroscopic observations and their reduction, including 
the treatment of host galaxy subtraction and flux cali- 
bration. In § |4] we consider our sample with respect to 
the broader set of SNe found by SNLS, ensuring it is a 
representative subset in terms of various observables, and 
discuss existing local UV spectra. In § [5l we undertake 
the detailed analysis. First we compare the UV spectra 
found in our sample with those found locally. We then 
examine the diversity of intermediate-redshift SNe la in 
various ways and correlate the UV variations with the 
light curves of the SNe and the properties of the host 
galaxies. We discuss these trends in terms of progenitor 
mechanisms in §[6] and examine the implications in terms 
of possible long term limitations of SNe la as probes of 
dark energy. We also present the mean phase-dependent 
SN la spectrum and its uncertainties for use in future 
work. 

2. SELECTION OF SNLS SNE lA 

Our SNe la are taken from the Supernova Legacy Sur- 
vey (SNLS), a "rolling" search for distant SNe with a pri- 
mary science goal of determining the avera ge equation- 
of-state parameter of dark energy, <w> (see lAstier et al.l 
2006). S NLS exploits th e square-degree Megacam cam- 
era (jBoulade et al |[2003D on the Canada- France-Hawaii 
Telescope (CFHT), and comprises repeat imaging in 
4 filters, g'r'i'z', of four deep 1x1 degree fields (see 
ISullivan et all 12006a' for the field coordinates), plus fur- 
ther u* imaging which is not timc-scquenced. Each field 
is imaged several (~ 5) times per lunation for 5-6 lu- 
nations per year. A description of the real-time search 
operations and the criteria for f ollowing SN ca i ididate s 
spectroscopically can be found in 'Sulliva n et all ()2006af) . 
Spectroscopic follow-up time for the essential work of ba- 
sic redshift measurement and SN type determination for 
cosmological analyses comes from major long-term pro- 
grams at the European Southern Observatory Very Large 
Telescopes (Basa et al., in prep; PI: Pain), the Keck ob- 
servatory (PI: P erlmutter) and the Gemini North and 
South telescopes (jHowell et al.ll2005l:IFronder et al.ll2007l 
PI: Hook). 

Observing time with the Keck-I telescope for the de- 
tailed spectroscopic study of individual SNe presented in 
this paper was scheduled from 2003 through 2005. The 
goal was to obtain substantially higher signal-to-noise 
ratio spectra than those required for SN typing (particu- 
larly in the rest-frame UV) , and hence we targeted SNe la 
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with a mean redshift lower than that of the SNLS as a 
whole (z=0.45 versus z=0.6) and used integration times 
~ 3-4 times longer for that redshift. 

To avoid spectroscopic screening of every candidate 
(results from the mainstream VLT and Gemini programs 
were typically not available at the time of the Keck ob- 
servations), a code was developed to constrain the red- 
shift, SN type and epoch from the available photometric 
data (see lSullivan et al1l2006al for details). In selecting 
candidates for this program, the following criteria were 
adopted: 

1. The SN photometric redshift lay at z<0.75, 

2. The predicted phase was such that the Keck ob- 
servations would be conducted prior to or close to 
maximum light. 

3. The light curve and color were consistent with a 
SN la. 

In practise, this last criteria is conservatively applied 
to minimize the risk of selecting against SNe la which 
may differ from the template used in our selection code, 
though this reduces the fraction of SNe la observed: By 
the end of 2005, after 20 nights of Keck time, 58 SN 
candidates were observed of which 36 were confirmed as 
SNe la (^^60%). A higher fraction of SNe la could have 
been obtained by a more rigorous application of the pre- 
selection technique; this argues that any bias introduced 
by our selection is small. 

Table [1] summarizes the key parameters for the SNe la 
sampled in the Keck campaign (an independent sur- 
vey of SNe IIP w as also undertaken during this period 
(jNugent et al.l200Q') .) We restrict these 36 SNe to a high- 
quality subset, removing 3 as spectroscopically peculiar 
(see ii l5.ip . 2 with a low-S/N Keck spectrum, and 5 with 
poor g' light curve coverage which precludes an accurate 
subtraction of the host galaxy from the SN spectrum 
(see Table [1]). These later 5 were all observed during the 
SNLS "pre-survey" when the light curve coverage in g' 
was less dense, and hence our photometrically calibrated 
SN spectrum will suffer a greater systematic uncertainty 
(see ij 13. 3p . Analysis of the remaining 26 SNe la forms 
the basis of this paper. Fig. [1] displays the redshift distri- 
bution of these high-quality SNe la and the breakdown 
in phase at which the Keck spectra were obtained. 

For comparison purposes, we will also make use of 
an additional sample of 160 spectroscopically-confirmed 
SNLS SNe la (without high quality Keck coverage) 
within the above redshift range for which good broad- 
band photometric data is available. This will be used as 
a control to demonstrate our SN la selection is unbiased 
relative to the larger SNLS population in terms of funda- 
mental properties such as light curve stretch, host galaxy 
type and location within each host. 

3. SPECTROSCOPIC DATA 

3.1. Observations 

All spectroscopic observations of SNLS SNe presented 
here were conducted with the 10 meter Keck-I telescope 
using the two -channel Low Re solution Imaging Spectro- 
graph (LRIS; [Okeiiin EUl). Although SNe la have 
broad absorption features, the need to monitor expan- 
sion velocities to reasonable precision as well as to lo- 
cate narrower diagnostic lines which would otherwise be 



blended in broader features defined the need for inter- 
mediate (600 1 mm~^) dispersion gratings. Access to the 
rest-frame UV over 0.2<z<0.7 was a prime goal given 
the sensitivity of features in this region to metallicity, 
temperature and photospheric expansion velocity. 

Taking these considerations into account, we adopted 
a spectroscopic set up using a 400 linemm"^ grating 
blazed at 8500A ( giving a dispersion of ~ 1.9 A pix ) on 
the red side and a 600 line grism blazed at 4000A (dis- 
persion ~ 0.6Apix~^) on the blue side. With a 5600A 
dichroic, this gave a contiguous wavelength coverage from 
the atmospheric limit on the blue side through to 9400A 
on the red side. On some of the earliest runs in 2003, 
we experimented with a 300 line grism blazed at 5000A 
(dispersion ~ 1.4Apix~^) on the blue side and a 600 
line grating blazed at 7500A (dispersion ~ 1.3Apix~^) 
on the red side with a 68OOA dichroic. This provided a 
lower resolution on the blue side, but a more optimal set 
up on the red side for the higher redshift SNe. However, 
we found the 68OOA dichroic did not permit accurate flux 
calibration on the blue side, because of residual contam- 
ination in the second order of diffraction. In very good 
seeing conditions, we used a 0.7" slit, otherwise a 1" slit 
was used. 

Ob servations were p erformed at the parallactic angle 
(e.g. 'Filip penkol 119821 ). Each individual exposure was 
typically between 1200-1800s, and the SN was dithered 
along the slit by 3-5" between exposures to aid with 
fringe removal in the data reduction. The telescope was 
refocused periodically during each night from which an 
estimate of the seeing was measured. Table [5] summarizes 
these spectroscopic observations. 

3.2. Data Reduction 

All spectroscopic data were processed using a pipeline 
developed by one of us (MS). This pipeline uses 
both standard Image Reduction and Analysis Facility 
(IRAF^") software as well as our own custom- written 
routines. LRIS is a two-arm spectrograph and different 
reduction techniques are required on each side as the red 
side suffers signiilcant fringing above 7000A. 

Our first step is to remove the overscan level on each 
amplifier of each CCD, and then subtract a master zero- 
frame constructed on the afternoon of each night of ob- 
serving to remove any bias pattern. Both sides are then 
divided by a normalized internal flat-field which removes 
pixel-to-pixel sensitivity variations. This division also 
corrects for the different gains of the LRIS CCD ampli- 
fiers. The 0.7" slit additionally requires division by a 
"slit-flat" or "illumination" flat-field as the slit illumina- 
tion is quite uneven, presumably due to milling defects 
in the slit. An additional complication is that the spatial 
position of the slit on the CCD can drift by several pix- 
els over the course of the night, so the slit-flat is shifted 
to best match each observat ion. Cosmic-rays a re then 
identifled using LACOSMIC ijvan Dokkum|[200l[) and re- 
moved via interpolation from neighboring pixels. 

We perf orm our sky su btraction following the tech- 
nique of iKelsonI (|2003D . which subtracts a two- 

^'^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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dimensional sky frame constructed from sub-pixel sam- 
pling of the background spectrum and a knowledge 
of the wavelength distortions as determined from two- 
dimensional arc comparison frames. For the blue side, 
we proceed directly to the spectral extraction after sky 
subtraction; on the red side we perform an additional 
fringe correction before extraction. 

To generate the master fringe frame, we first divide 
every sky-subtracted frame by the two-dimensional sky 
frame that was subtracted and that caused the fring- 
ing, resulting in a two-dimensional map of the fringe 
strengths per photon for each exposure. All individual 
fringe maps taken in a given configuration are averaged 
with sigma-clipping to form a master fringe frame, mask- 
ing out regions in each exposure that contain object flux, 
and weighting each frame by the exposure time of the ob- 
servation. As the observations are dithered, the objects 
are located on a different part of the CCD in each frame, 
and the fringe strength at every pixel is therefore sam- 
pled across an entire night. Finally, we boxcar smooth 
the master fringe frame to increase the signal-to-noise 
ratio in the fringes, which typically extend over several 
pixels. This smoothed, master fringe frame is then scaled 
and subtracted from each red-side exposure at the stage 
prior to sky-subtraction, and the sky-subtraction then 
repeated on the fringe-subtracted images. 

The two-dimensional frames are transformed to a con- 
stant dispersion using comparison arc lamp images. Be- 
cause the objects under study typically have a reason- 
able signal-to-noise ratio, the spectral extraction was per- 
formed by tracing the object position on the CCD and 
using a variance- weighted extraction in a seeing matched 
aperture, equal to 1.25 times the seeing. An error spec- 
trum from the statistics of the photon noise is also ex- 
tracted at the same time. The wavelength calibration of 
each extracted spectrum is then tweaked using the po- 
sition of the night-sky lines to account for any drift in 
the wavelength solution during a given night. We then 
perform an approximate telluric and relative flux calibra- 
tion using spe ctrophotometric standard stars taken from 
iMassev et al.l ((1988), and correcting for typical Mauna 
Kea atmospheric extinction using the extinction law of 
iKrisciunas et al.l (|1987t l and the effective airmass of each 
exposure. We note that flux calibration errors introduced 
during this process that vary smoothly as a function of 
wavelength will be corrected by our calibration to the SN 
lightcurves described in the next section. 

Following optimal combination of the different expo- 
sures and relative flux calibration of the blue and red 
components, our final step is to scale the two sides into 
a single calibrated spectrum. We first match the flux 
across the dichroic by defining narrow box filters either 
side of the dichroic and scale each spectrum appropri- 
ately. We then use a weighted mean to combine the 
spectra, in the process re-binning to a constant 2A per 
pixel resolution. The result is a contiguous science spec- 
trum, together with an error spectrum representing the 
statistical uncertainties in the flux in each binned pixel. 

Where possible, we measure redshifts from lines in the 
host galaxy spectrum. In some cases the host galaxy was 
too faint, or the SN too isolated from the host galaxy, to 
measure a redshift from the galaxy lines, and the redshift 
was determined from t emplate fits to the SN spectrum 
(e.g. lHowell et al]l2002[ ). As a final step prior to scientific 



analysis we interpolate across the positions of host galaxy 
emission features (e.g.. Ha, H/?, O II, etc.) to remove 
these from the spectra. 

The SNe la themselves were identified by a lack of hy- 
drogen in their spectra, combined with broad (several 
thousand kms^^) P Cygni lines of elements such as Si II, 
S II, Ca II, Mg II, and blends of Fe peak lines. The clas- 
sification of SNe is discussed extensively i n the litera- 
ture; seejFilippcnko ( 1997) for a. review, and I Hook et al.l 
i2005), Li dman et all (120051). iMatheson et al.l (120051 . 
and Howell et (|2005[ ) for issues associated with high- 
redshift SN classification. SNLS as a who le adopts the 
numerical classification scheme presented in lHowell et al.l 
(2005) for typing SNe: a confidence index (CI) that a 
given spectrum is of a SN la. Three classifications denote 
SNe la: "Certain SN la" (CI:5; type denoted as SN la), 
"Highly probable SN la" (CI:4; type denoted as SN la) 
and "Probable SN la" (CI:3; type denoted as SN la*). In 
all of the spectra in this paper, the relatively low-redshift 
compared to the SNLS sample as a whole, coupled with 
the long integration times, results in high signal-to-noise 
ratios ensuring the classifications are unambiguous. All 
SNe la presented here are CI=5 on the SNLS scheme. 

3.3. Spectrophotometric Calibration and Host Galaxy 
Subtraction 

Although standard spectrophotometric stars were ob- 
served periodically on each observing night, these provide 
only an approximation to a flux calibrated spectrum. 
Complexities such as differential slit-losses (as the see- 
ing is wavelength dependent) and uncertain atmospheric 
extinction corrections require a more rigorous treatment. 
Furthermore, a key issue in the interpretation of SN spec- 
tra is the removal of underlying contaminating flux from 
the host galaxy, which for some high-redshift spectra can 
be significant. We adopted a procedure for removing the 
host galaxy spectrum that takes advantage of the pho- 
tometric u*g'r'i'z' host-galaxy data available from the 
SNLS, as well as the near-simultaneous SNLS g'r'i'z' SN 
light-curve photometry, ensuring that each spectrum is 
carefully corrected to be consistent with the appropriate 
broad-band SN colors. 

We begin by measuring the surface brightness of each 
galaxy (in flux units per square arcsecond) at the posi- 
tion of the SN in a series of small apertures (with radii of 
1-8 pixels) for each of the u*g'r'i' z' filters. These fluxes 
were measured from deep stacks with no SN light present, 
constructed for the pur pose of studying the p roperties of 
the SN host galaxies (jSullivan et all l2006bf) . The im- 
age quality of these stacks was around 0.8" in each fil- 
ter, similar to the seeing for the majority of our LRIS 
spectral observations. The amount of host galaxy light 
present in each SN spectrum is then estimated by inter- 
polating the galaxy surface brightnesses at the seeing of 
the observation, and then multiplying by the product of 
the slitwidth (5*^^) and the spectral extraction aperture 
(1.25xseeing; § 13 2p . giving a series of u*g'r'i' z' fluxes 
representing the host galaxy contamination in our spec- 
tra. We then fit a series of smooth galaxy spectral energy 
distribution (SED) templates to these flux data to esti- 
mate the contaminating host galaxy spectrum. These 
templates are generated by the galaxy spectral synthesis 
code PEGASE.2 (|Fioc fc Rocca-Volmerangelll997l . ri999.) 
rebinned to the resolution of our SN spectra. We inter- 
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polate between neighboring galaxy spectra to find the 
best-fitting template. The resulting spectrum represents 
the estimate of the contaminating host galaxy continuum 
in every SN spectrum. 

We next use the g'r'i' z' SNLS light curves and their re- 
sulting light-curve parameterization (using the lightcurve 
fitter SiFTO; Conley et al., in prep) to interpolate the SN 
flux at the time of each spectral observation from nearby 
light-curve observations. Due to the rolling-search na- 
ture of the SNLS and the dense light curve coverage, 
most of the spectra in our main sample have g' photo- 
metric data within 2 rest-frame days of the spectroscopic 
observation, and all have data within 4 days (Table [1]), 
hence the interpolation required is small. 

We next estimate the amount of SN flux that went 
through the LRIS slit and was inside the spectral ex- 
traction window of 1.25 x seeing, and hence is present in 
the reduced spectrum. We model the SN PSF as a two- 
dimensional Gaussian with a = seeing/2.355 and inte- 
grate this Gaussian through an rectangular aperture of 
5^1/ X (seeingx 1.25), where a is calculated at the effec- 
tive wavelengths of the g'r'i' z' filters adjusting the see- 
ing as A~°'^. This estimates the fraction of the total SN 
flux that passed through the slit, mimicking the inclu- 
sion of differential slit-losses. These adjusted SN fluxes 
are added to the host galaxy fluxes to generate combined 
"SN-|-host" g'r'i' z' fluxes, the estimate of the amount of 
flux that passed through the slit. 

We correct our combined, contiguous and flux- 
calibrated spectrum to have the same colors and absolute 
flux level as dictated by these "SN-|-host" fluxes using a 
smooth interpo lating multiplicative spline function (see 
iHsiao et aLllBoTl . The size of the correction made to the 
observed spectrum is usually <10% and is invariably a 
monotonic function of wavelength. We then subtract our 
best-fitting host galaxy template, and adjust the result- 
ing spectrum to have the correct colors of the SN on the 
night of observation (again using a spline function) , giv- 
ing a SN spectrum with the correct relative and absolute 
flux calibration. This flnal multiplicative adjustment is 
small and essentially corrects the subtracted spectrum 
for differential slit losses. We note that this technique 
would not be possible without either the u*g'r'i'z' host 
galaxy fluxes, or the densely-sampled g'r'i' z' SN la hght 
curves. 

We estimate the uncertainty in the host subtraction 
process using a Monte-Carlo simulation. We repeat the 
subtraction for each SN 500 times, but adjusting the dif- 
ferent variables (seeing, host-galaxy fluxes, SN fluxes) ac- 
cording to their uncertainties and assuming normal dis- 
tributions. We assume a 15% uncertainty in the seeing, 
and also simulate errors in the centering of the SN in the 
slit, introduced during target acquisition or by telescope 
tracking errors, of ±0.15" {la). Where appropriate, we 
use these 500 simulated spectra to estimate errors due 
to host galaxy subtraction in derived quantities (such as 
spectral colors) in later sections. We also carry forward 
in our analysis the statistical error arising from the pho- 
ton statistics in our spectra 13. 2p . 

Our final task is to correct the SN la spectrum 
for the effects of extinction in the Milky- Way. To 

correct for Galactic exti nction, we use the dust- 

maps of .Schlegcl, Finkbeiner. fc Davis! (|1998D and a 
ICardeUi. Clavton. fc MathisI (|1989l ) extinction law. We 



defer the more complex discussion of extinction in the 
SN host galaxy and intrinsic SN color variations until 
later. 

We illustrate the results of this process for all 26 of our 
high quahty SNe la in FigM^ 

4. METHODOLOGY 

We now turn to the two key questions our dataset is 
designed to address: 

1. Has there been significant evolution in the mean 
SN la spectrum since z ~0.5 - a period of 5 Gyr 
corresponding to the past 37% of cosmic history? 

2. To what extent are there intrinsic spectral vari- 
ations in our sample? Do these give indications 
that SNe la are a more complex population than 
assumed in cosmological studies, and might these 
variations limit the use of SNe la in future cosmol- 
ogy experiments? 

4.1. Theoretical Predictions 

We will begin to address these questions by exploring 
the role that our rest-frame UV data offers as a proxy 
for progenitor metallicity. The possible effects of progen- 
itor metallicity, either as a function of redshift through 
galaxy evolution, or as a function of environment or host- 
galaxy luminosity, can be explored through blanketing 
and wavelength dependent features in the rest-frame UV 
corresponding to AA 2900-3500A. At z>0.2, this portion 
of the rest-frame spectrum is well-sampled in our LRIS 
data. In making the assumption that our UV data are 
indicative of metallicity effects, it should be realized we 
are relying largely on theoretical studies; the UV spec- 
tral region is poorly explored observationally. First we 
r eview the various theoretical expectation s. 

iHoflich. Wheeler, fc Thielemann! () 19981 ) argue that di- 
rect traces of the progenitor metallicity can best be seen 
in the unburned SN layers which are only observable sig- 
nificantly before maximum light. However, they also pre- 
dict that an increase in progenitor metallicity will cause 
an increase in the amount of ^^Fe synthesized in the ex- 
plosion, and this will result in an increase in line opacity 
in the UV region which may be observable at maximum 
light. The degree of mixing in the explosion complicates 
the interpretation, however. The net effect predicted by 
Hoflich et al. is that an increased metallicity will result 
in an increase in the UV pseudo-continuum at maximum 
li ght. 

Lentz et aLl (|2000r) start with the flndings of 
IHoflich et all (|199aD and examine the spectroscopic im- 
plications in greater detail. Using the starting isotopic 
distribution of the W7 model (Nomoto et al. 1984), they 
simultaneously change the progenitor metallicity in the 
unburned C-l-0 region and increase the amount of ^^Fe 
in the partially burned region. They find two effects - a 
shift in the wavelength of UV features redward with de- 
creasing metallicity, and a simultaneous increase in the 
level of the UV pseudo-continuum. They argue that as 
metallicity decreases so the line opacity decreases with 
the result that lines form deeper in the atmosphere, and 

The online version contains all spectra; the printed version 
contains two examples. 
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therefore from a lower velocity region. A lower metal- 
licity also decreases the cooling producing a higher tem- 
perature and bluer colors. We note that such a color 
difference is opposite to the effect predicted by Hoflich et 
al. However, Lentz et al. caution that the overall UV 
flux level is not necessarily a good indicator of metal- 
licity, as it is dependent on many variables such as the 
t emperature, density and velo c ity of the C+O layer. 

iTimmes. Brown, fc TruranI (|2003D argue for much 
more dramatic changes in SN la physics with metallicity. 
They argue that synthesized ^^Ni mass should be linearly 
proportional to progenitor metallicity. Since the decay of 
^^Ni drives the luminosity, SNe la in high metallicity en- 
vironments should be less luminous. This is because stars 
from higher metallicity environments will end up with 
larger mass fractions of ^^Ne and ^^Fe after helium burn- 
ing. Since these isotopes have excess neutrons, the au- 
thors argue that in these cases fewer radioactive elements 
are produced during the process of burning to nuclear 
statistical equilibrium during a SN la. Note that these 
results are in sharp contrast to other studies which found 
no si gnificant increase in^^Ni with increasing metallicity 
(e.g. iHofiich et aI1ll998t llwamoto et"allll999f ) 

Finally, we note that various authors predict that the 
SN la rate should be affected by metallicity, although 
they do not make explicit prediction s about the re- 
sultin g effects on SN la properties. iKobavashi et al.l 
(|1998D argue that in very low metallicity environments 
{[Fe/H] < — 1) the white dwarf wind that they beheve 
is essential for producing SNe la will b e inhibited, thus 
leading to fewer SNe la. Lang er et al.l (^000) find that 
metallicity differences should alter the range of progeni- 
tor masses that produce SNe la. 

In summary therefore, theory cannot yet offer us a 
clear consensus as to the effects of metallicity on the 
UV properties of SNe la. Indeed, there is disagreement 
not only about which effects are the most important, but 
also about the sign of any possible effect. This is a very 
challenging theoretical problem hindered by correlations 
between the wanted effect of metallicity and other cor- 
relations such as the viability of certain progenitor sys- 
tems, the explosion mechanism and radiative transfer in 
an atmosphere under a variety of mixing conditions. Full 
simulations of all these effects may soon become feasible, 
but substantial campaigns will still be needed to track 
statistical shifts with metallicity. We therefore under- 
take an empirically- motivated study of the properties of 
supernovae in the UV, and note where the observations 
agree or disagree with certain theoretical studies. In this 
respect, it is helpful to have some theoretical basis for 
making the empirical measures. 

Fig. Bl illustrates, for th e range of metallicities sam- 
pled bv lLentz et al.l (|2000( ). that the bulk of the spectral 
changes would clearly be detectable within our rest-frame 
wavelength range. The blanketing and wavelength shifts 
apparent in Fig. [3] suggest the use of both photometric 
and spectroscopic diagnostics. As high quality UV spec- 
tra are in scarce supply at low redshift, we can best make 
meaningful comparisons using rest-frame U — B measure- 
ments. As our spectra do not always extend to the full 
redward end of the standard B filter, we will also adopt 
a top-hat filter extending from 4000 to 4800A which we 
will refer to as b. Pseudo-photometric ultraviolet mea- 
sures can be constructed from our spectral data using two 



top-hat filters, UVl and UV2, centered around regions 
where the strongest trends are predicted. Although these 
various photometric measures degrade the information 
content of our spectra, they have the considerable benefit 
of being consistent over a wide range in redshift without 
making any assumptions about the k-correction. The 
metallicity-induced spectral shifts discussed above can 
likewise be tracked by measuring the rest wavelengths of 
two UV features at AA2920 and 3180 A . We wiU refer 
to these indicators as Ai and A2. Both the photometric 
and wavelength features introduced above are marked in 
Fig. El 

4.2. Suitability of the Keck SNe la Sub-sample for 
Detailed Studies 

If our conclusions concerning the properties of z ~0.5 
SNe la are to be relevant for the role of SNe la as cosmo- 
logical probes, our sample of SNe la should be a represen- 
tative subset of those being used for probing dark energy. 
Accordingly, first we will verify whether this is the case 
by comparing the distribution of properties with those 
for the parent SNLS SN la sample. This parent sampl e 
(N=160) includes those presented bv lAstier et al.l (|2006D . 
as well as further events up to May 2006, restricted to 
lie within our chosen redshift range 0.15< z <0.7. This 
is the redshift range over which SNLS suffers least from 
Malmquist-type effects in the properties of SNe chosen 
for spectroscopic follow-up, and the broad SNLS popula- 
tion is reasonably unbiased oyer th is redshift range (e.g. 
lAstier et al.ll2006l: iNeill eT^\20M) . 

Fig. |4] compares the distribution of SN la stretch, lo- 
cation and host galaxy properties for the two samples. 
There is no evidence that our selection criteria (§2) have 
significantly biased the SN la population chosen for in- 
tensive study relative to the larger SNLS population. 
Kolmogorov-Smirnov (K-S) tests indicate the Keck and 
SNLS sample distributions in all four parameters are con- 
sistent with being drawn from the same distribution; the 
exception is the stretch distribution where the K-S test 
indicates the distributions are possibly different due to 
a marginal deficiency of low stretch events in the Keck 
sample. This can be easily understood as these events 
are both fainter (making the Keck spectroscopic obser- 
vations more demanding) and have faster light curves 
(making a pre-maximum or at maximum light observa- 
tion less likely, a criterion when selecting SNe for Keck 
foUowup; §12) than higher-stretch events. Otherwise, the 
larger scale environs in which our SNe la occurred seem 
quite representative of the parent population, for exam- 
ple in the specific star formation rate, the separation from 
the host galaxy, and (not shown in Fig. |3|) the absolute 
magnitude of the host galaxy. 

4.3. Local UV Spectra 

In order to examine possible evolutionary effects we re- 
quire a local b aseline UV spectrum. Fo r this, we turn to 
the analysis of ' Nugent. Kim, fc Perlmutteri (j20l)3 ) which 
provides a phase-dependent template SN la spectrum 
which has become the standard reference used exten- 
sively, for example, in the calculation of cross-color k- 
corrections in deriving SNe la luminosities at high red- 
shift. Further details of its construction can be found in 
the original article. Note t hat though improved spectral 
templates now exist (e.g. iGuv et all 120071 : iHsiao et"al] 
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l2007f) , these typicaUy exploit high-z spectra in their con- 
struction, whereas Nugent et al. use solely local SN la 
data. 

The UV portion of this mean local spectrum is based 
on very few local ev ents studied b y eithe r the lUE or 
HST (see Table 1 in iNugent et~all (|2002f )). and signif- 
icant uncertainties remain at the shortest wavelengths. 
Unfortunately, unless STIS is resuscitated on the up- 
coming HST Servicing Mission, there is no immediate 
prospect of improving this situation. A particular con- 
cern is the fact that within a week of peak brightness 
only three SN la (comprising of a total of five spectra) 
contribute to the UV template. They are SNe 1981B, 
1990N and 1992A. These SNe have stretches of 0.929, 
1.074 and 0.819 measured using SiFTO, i.e. within the 
range sampled here (Fig. 2]). However, their hosts, classi- 
fied as SAB(rs)bc (NGC 4536), SAB(rs)bc (NGC 4639) 
and SAO (NGC 1380) respectively, are possibly biased 
towards large early type systems. 

5. RESULTS 

5.1. Tests for Progenitor Evolution: the Mean UV 
SN la Spectrum 

We begin by comparing the mean spectrum of the 
int ermediate-redshift K eck SN la sample with the lo- 
cal iNugenFetlai] (|2002t l template. Earlier work compar- 
ing the properties of low and high-redshift S N la spec- 
tra fe .g.lCoil et al.ir2 000: Hook et al. 2005; Bl ondin et all 
120061 iGaravini et al.ll2007l : iBronder et al.ll2007^ has con- 
centrated largely on comparisons of optical spectral fea- 
tures involving lower-S/N spectra. In general, these stud- 
ies find that equivalent widths and ejection velocities of 
optical spectroscopic features measured on high-redshift 
SNe la have similar distributions to those observed lo- 
cally. No study has yet found evidence for evolution in 
the optical spectroscopic properties of SNe la with red- 
shift. 

In this paper, we concentrate on the properties of the 
mean UV spectra. As discussed earlier, models sug- 
gest metallicity variations will have greatest effect in this 
wavelength region and thus our sensitivity to evolution- 
ary effects should be maximized. We separate the distant 
SNLS events into three groups (Fig. [1]): those observed 
prior to maximum light [t <-4 rest-frame days), those 
observed close to maximum light (-4< t <-|-4 days), and 
those observed after maximum light (t >-|-4 days). We 
exclude 3 spectra, labeled (f) in Table [U which we clas- 
sify as peculiar. One is the super- Chandrasekhar event 
SNLS-03D3bb ()Howell et al.l I2006D and the other two 
(SNLS-05Dlhk and SNLS-03D4cj) ha ve peculiar spec- 
tra s imilar to that of SN1991T-type (Filip penko et al.) 
Il992f) . (A further 7 events, labeled (d-e) in Tabled were 
discarded as described in § [2l) At this stage, our com- 
parisons will involve spectra uncorrected for host galaxy 
extinction and/or intrinsic color variations. 

Our 26 high-quality SNe la consist of 8 pre-max, 15 
maximum light, and 3 post-max spectra. In creating a 
composite spectrum, flux-calibrated spectra are normal- 
ized using the top-hat b filter introduced earlier which 
runs from AA4000-4800A (see Fig. El). We use a clipped 
mean and determine the error on the mean by bootstrap 
resampling. To facilitate a meaningful comparison we 
color-adjust the local spectrum so that its U — B color 



matches that of the mean high redshift sample. Figs. \5\ 
and [6] show the maximum- light and pre-maximum high- 
z mean spectra compared against the local template at 
similar phases. For convenience, we present in Tables [3] 
andlH the early and maximum light spectral energy dis- 
tributions with upper and lower 90% confidence limits. 
For the local template, due to the very small number 
of SNe la with rest-frame UV spectra (§ I4.3p . a similar 
bootstrap resampling technique to asses the uncertainty 
in the local template is not meaningful. 

Examining Fig. [5l it is clear that the spectrum at 
z ~0.5 is reasonably similar to its local equivalent. The 
agreement, particularly longward of4000A, IS reassurmg, 
indicating that there has been no significant evolution in 
the mean SN la spectroscopic properties over the past 
5 Gyr. However, some systematic departures are seen, 
for example, in the Si II/Co II doublets at ~4150A and 
in the positions of the UV features around 2900-3200A. 

Do these UV changes represent systematic differences 
in the progenitor properties between z ~0.5 and 0? In 
this respect it is helpful to consider the error in the z ~0.5 
mean spectrum (shaded in Fig. O. Here we see that the 
local spectrum lies within the 90% confidence boundary 
in the distribution of z ~0.5 spectra. Thus it is quite pos- 
sible that the UV differences seen between high and low 
redshift arise as a result of statistical variations occurring 
within a non-evolving population, especially considering 
the number of local SNe la with rest-frame UV coverage 
is very small. 

Turning to the pre-maximum light spectra (F ig. [6]), 
both Hofl ich et al.l (^1998) and Lentz et all (|2000D argue 
that differences in the UV spectrum due to changes in the 
progenitor metallicity should be larger at earlier times, 
where the unburned outer layers of the white dwarf play 
a larger role in shaping the appearance of the spectrum. 
Although the agreement does indeed seem to be poorer, 
the local template is considerably uncertain at this phase; 
we conclude that the most likely explanation for the dif- 
ferences we see is that the local template is unrepresen- 
tative in the far UV. 

However, the differences that we do see in the UV data, 
whether due to evolution or simply as a result of intrin- 
sic scatter in the p opulation, are mu ch larger than the 
effects predicted bv lLentz et al.l ()2000l ) for quite substan- 
tial changes in progenitor metallicity. It should be noted 
that in these models, luminosity and kinetic energy were 
fixed, whereas our observations presumably span a con- 
siderable range in both of these variables. Lentz et al. 
also predict systematic wavelength shifts of the diagnos- 
tic features identified as Ai and A2 in Fig. [H the study 
of which should be independent of reddening. We will 
defer consideration of these spectroscopic features until 
§[5X1 

At longer wavelengths, the local data is more reliable. 
Here, several features of the pre-maximum spectrum de- 
serve comment (Fig. [6]). Compared to the low-z template, 
the high-z average spectrum has a shallower S II 5400A 
feature, and a shallower Si II 4130A feature. In addition 
the 37OOA Ca-^Si feature is narrower and shows a greater 
degree of splitting into the separate Ca and Si compo- 
nents in the high-z c omposite spectrum . More high-z 
SNe la belong to what lBranch et aTl (|2006( ) call the "shal- 
low Si" group — SNe such as SN 1991T, which have a 
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narrower zone of intermediate mass elements. These re- 
sults agree with the findings of iLi et all (|2001l ) who note 
that, while a significant fraction of SNe la show such be- 
havior at early times, they become more n ormal at max- 
imum light. This supports the claims of I Sullivan et all 
p006b) : at high redshift, where the fraction of star form- 
ing hosts is higher, more "prompt" SNe are found and 
these tend to have broader light curves. 

The lack of a comprehensive local sample with compa- 
rable quality to the Keck data extending into the UV 
clearly limits any precise tests for spectral evolution. 
This impasse raises the question of whether it might be 
more profitable to search for evolutionary trends within 
the Keck sample itself. Fig. [7] shows the mean maximum 
light spectra compiled for two Keck sub-samples, split 
into two redshift ranges at z=0.50. The mean redshifts 
of the two subsamples are 0=0.36 {N=7) and 2=0.58 
(N=8). The agreement is striking, with only marginal ev- 
idence for differences in blanketing at the shortest wave- 
lengths. Further progress in this respect may be possi- 
ble with the ongoing SDSS SN la survey (which samples 
z ~0.2) and the higher-redshift SNLS SN la spectra at 
z ~ 0.6-1.0. 

In summary, no convincing evidence is available for 
a systematic change with redshift in the UV spectra of 
SNe la. The bulk of the differences seen probably arise 
from a natural dispersion within the sample and it seems 
reasonable to assume this dispersion is present at all 
epochs. In the next section we explore the possible origin 
of this UV dispersion. 

5.2. Reality of the UV Dispersion 

In this section we first consider the physical reality 
of the dispersion seen in the UV spectra at maximum 
light. The bootstrap-resampled variations (Figs. [5] and 
[5]) indicate a greater degree of scatter at A < 3700A than 
in the optical around A ~ 5000A. However, both the 
effects of dust extinction and the size of the calibration 
uncertainties are unaccounted for in these comparisons. 
Could the increased dispersion simply arise as a result 
of one or both of these effects and not represent a true 
physical variation amongst our SNe? 

We first discuss the effects of dust extinction. 

5.2.1. Effects of Dust Extinction 

An explanation for the UV variations is the pres- 
ence of varying amount of interstellar dust, either in 
each host galaxy or along the line of sight to the ob- 
server. Correcting the properties of SNe la for extinc- 
tion is currently one of the most pernicious problems 
in the ir use as cosmolog ical probes (e.g., see discus- 
sion in IConley et al.ll2007ri. Earher a n alvses (iRiess et al.l 
19961 119981: iPerlmutter et all 119991: iRiess et all 120001 : 
Sullivan et al.ll2003HKnop et al.b oOSh. based on the first 
generation of SN la surveys, examined the range of likely 
extinction using both SN colors and host-galaxy mor- 
phologies as markers. Generally only modest levels of 
extinction {Ay < 0™.2) were seen at high-redshift across 
the full range of host types, with no evidence of any 
systematic change with redshift in the samples used for 
cosmological analyses. 

In the case of the present SNLS sample, the key di- 
agnostic is the rest-frame B — V color estimator "c" 



(|Guv et al.ll2005l .l2007f). essentially theB-V color mea- 
sured at maximum light. We wish to test whether dust is 
the primary cause of the UV variations seen in the Keck 
spectra. This can be done by considering, statistically, 
correlations between the UV spectra and c. Redden- 
ing measures for each SN are not directly output by the 
SNLS analysis since c is treated as an empi rical variable 
when performing the cosmological analysis (|Astier et al.l 
l2006( l: no physical model is assumed for the color varia- 
tions. However, we can readily test the hypothesis that 
the bulk of the UV dispersion arises from reddening that 
is linked to the B — V variations. 

We investigate two methods for color-correcting our 
SNe la spectra. The first is to use a standard reddening 
law (ICardelli e t al.l ()1989[ ). CCM) and extinction-correct 
each individual Keck spectrum using that SN's c mea- 
surement and a value oi Rb =4.1, suitable for standard 
Milky- Way type dust. The results from this process are 
not particularly sensitive to either the exact B — V color 
used as the fiducial/zeropoint SN la color, nor to the ex- 
act value of Rb used. Our secon d approach is similar but 
uses the SALT (jGuv et al.ll2005[) color law together with 
the c estimated from the light curve fits. As this color law 
was trained on local SN spectra and colors, it includes 
not only the effect of dust extinction but may also par- 
tially account for intrinsic variations between SN optical 
colors and UV fluxes unrelated to dust. The CCM and 
SALT laws have a similar form in the optical, but differ in 
the near-UV, where the SALT law imp lies more "extinc - 
tion" in the UV than the CCM law ()Guv et al.l[200l . 
To test the validity of these approaches we split those 
Keck spectra studied at maximum light into two sub- 
samples according to whether the rest-frame maximum 
light optical color is redder or bluer than a c=-0.057. 

Examining the mean UV spectra for the two subsets 
(Fig. [8|), we find that, before any correction, optically 
bluer SNe la do indeed have more flux in the UV as 
expected if they are less extincted. When the spectra 
are individually color-corrected using a CCM law and 
then combined by subset, the mean UV spectra agree 
better though the entire color difference is not corrected 
shortward of 3500A. The approach using the SALT color 
law performs better in this UV region, and provides a 
good agreement across our entire wavelength range, but 
is unable to correct all the differences seen. Given the 
slightly superior performance of the SALT color law, we 
adopt this approach for color-correcting our spectra in 
the remainder of this paper. 

5.2.2. Effects of Calibration Uncertainties 

Our host galaxy subtraction and spectral calibration 
procedures are described in § 13.31 For each host- 
subtracted and photometrically-calibrated SN spectrum 
we also record a series of 500 Monte-Carlo simulated 
spectra generated by randomly varying the SN fluxes, 
host galaxy fluxes, seeing and SN position in the slit ac- 
cording to the observational errors. 

For each of the SN studied, we estimate the 1-cr un- 
certainty as follows. In each wavelength bin for a given 
SN, we order the 500 random spectra and find the range 
in flux that encompasses 68% of the population and as- 
sign this flux range as the l-tr error in this wavelength 
bin. This is repeated in every wavelength bin for all 
the spectra, providing the error from the observational 
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uncertainties for each SN spectrum as a function of wave- 
length. 

We can then assess whether the scatter that we see 
about the mean spectrum is significant. For every spec- 
trum, we calculate the residual from the mean spectrum 
and divide by the appropriate error spectrum. The re- 
sult is the deviation from the mean for each spectrum 
in units of sigma. Smoothed versions of these are plot- 
ted in Fig. [21 the top panel shows the deviation for the 
uncorrected spectra, the middle and lower panels show 
the deviations for spectra corrected using the CCM and 
SALT laws as described in § 15.2.11 

In general, the scatter is markedly larger in the UV 
(A < 4000 A) than in the optical region. In the opti- 
cal, the spectra appear well-behaved, with most of the 
deviations <2(t, consistent with a mean spectrum show- 
ing little intrinsic dispersion over the broad wavelength 
range. 

These tests and those of the previous section demon- 
strate that most of the difference in the average contin- 
uum levels seen between the mean UV spectra of opti- 
cally blue and red SNe la can be corrected using exist- 
ing techniques such as the CCM dust extinction law or 
the SALT SN la color law. However, these techniques 
can only reduce, and not eliminate, the scatter observed, 
which is significantly larger than that expected from a 
consideration of the observational uncertainties. We thus 
conclude that the bulk of the dispersion is intrinsic to our 
SNe and investigate its origin in the next sections. 

5.3. Physical Correlations within the UV Dispersion 
5.3.1. Photometric Comparisons 

Given we have established that the UV variations are 
larger than those expected from photometric errors and 
are still present after attempts to correct for color dif- 
ferences, it is important to consider the extent to which 
the intrinsic dispersion we see is consistent or otherwise 
with that observed locally. A lthough no comp arable local 
spectroscopic dataset exists, iJha et all (|2006[ ) have ana- 
lyzed a homogeneous photometric database of 44 local 
SNe la and examined the U — B color dispersion at max- 
imum light in terms of both the light curve stretch factor 
s and redder colors such as B ~V . As in the higher red- 
shift data, there is still the complication of separating the 
effects of host galaxy dust extinction and intrinsic color 
effects. In the Jha et al. study, host galaxy extinction 
estimates were available for a subset of the data yielding 
intrinsic U — B and B — V colors. The increased intrinsic 
dispersion seen in the U — B vs. stretch relationship with 
respect to that m. B — V can thus be attributed primar- 
ily to an increased dispersion in the U band of ajj 2±0.12 
(c ^l2%). 

lAstier et al.l ()2006[ ) also address the question of the 
intrinsic dispersion in U by comparing the photometric 
properties of that subset of SNLS SNe la for which some 
of the g'r'i' z' bands map conveniently onto rest-frame 
U BV . They estimate the quantity A which is the dif- 
ference between the U band flux at maximum light pre- 
dicted from a suitably-chosen triplet of observed filters 
drawn from g'r'i' z' and the actual observed U magni- 
tude. As the machinery adopted by Asticr et al. includes 
provision for a color-stretch relation, this is not quite the 
same test as that adopted by Jha et al., although the 



local scatter found in At/a is comparable (0.12). 

We compare our intermediate redshift U — 6-stretch re- 
lation with that derived from the local sample in Fig. [TUl 
In this figure all colors are color-corrected using a SALT 
color law and the method of § 15.2.11 To derive the rest- 
frame U — b and c measures for the low-z sample, we re- 
perform the light curve fits to the SN photometry using 
the same fitter used for the SNLS sample (SiFTO; Conley 
et al. in prep). Clearly, at maximum light the trends are 
very similar. As discussed in § [51 the Keck sample (and 
the SNLS sample in general) is slightly under-represented 
in low-stretch SNe la compared to the local data (Fig. [4]), 
so the color-stretch relation found in previous studies is 
less evident. Crucially, there is no evidence that the col- 
ors derived directly from the Keck spectra show a larger 
photometric dispersion. For events with stretch s >0.85 
we find ([/ — 6)=-0.452 with a standard deviation of 0.132 
for 15 Keck SNe la c.f. -0.436 (0.167) for 41 local events. 
The small difference in mean color probably arises from 
the different selection criteria, as Jha et al. are able to 
include events with greater extinction. Note there are 
small differences in the definition of "maximum light" ; 
the Keck spectra range over -4 to -1-4 days which will 
add an additional scatter to this sample in the presence 
of color evolution with phase. 

A similar behavior is seen in the {UV2 — 6) - stretch re- 
lation; a modest dispersion together with a weak correla- 
tion with stretch and a slight phase dependence (the def- 
initions of UVl and UV2 are shown in Fig[3]). However, 
a marked increase in scatter is seen in the {UVl — b)- 
stretch relation (as expected if the UVl filter is the 
more sensitive to variations in progenitor composition) 
together with a stronger evolution with phase. Although 
the UVl filter is the shortest wavelength probed by the 
Keck spectra, this dispersion is considerable even within 
a phase bin, and is unlikely to arise from spectropho- 
tometric errors given the tight dispersion observed in 
the UV2 — b relation. Fig. [TUl also shows weak corre- 
lations between color and stretch in the sense that the 
high stretch (brighter) events are bluer tha n those of 
low s t retch (fainter) (e.g. Phillips ct al. 1999; Knop et al.l 
[2003t iGarnavich et al.ll2004l) . 

Clearly as the SNe la develop they display an in- 
creased dispersion in the wavelength region where metal- 
dependent features are expected to be most prominent. 
This may indicate that the metallicity of the outer, un- 
burned layers, seen only at early times, is n ot as impor- 
tant a s expected bv lHoflich et all ()1998f ) and lLentz et al.l 
(|2000D . Instead line blanketing from iron-peak elements 
may play the dominant role in affecting the UV opac- 
ity. This may affect maximum-light and post-maximum 
spectra more than for pre-maximum spectra. At early 
times the photosphere has not yet receded into the bulk 
of the iron-peak elements, and, as time progresses, ^^Ni 
synthesized in the explosion decays into ^^Co and subse- 
quently into ^^Fe. These additional isotopes may provide 
the increased line blanketing seen in the spectra. 

5.3.2. Spectroscopic Features 

Next we co nsider trends in the spectroscopic features 
introduced bv lLentz et al.l (|2000i ) as possible metallicity 
diagnostics; these are marked as Ai and A2 in Fig. [3l 
The wavelengths of these features have the advantage of 
being independent of any dust correction. If progeni- 
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tor evolution is largely driven by metallicity effects, we 
might expect significant shifts to occur in these features 
between the local and high-redshift data. 

In fact, a cursory examination of our mean spectra 
(Figs. [5] and [6]) shows that the features do not appear 
to shift significantly in the mean. However this could be 
due to an inadequate or unrepresentative local data. We 
have seen that this is a significant limitation. 

Accordingly, we measured the wavelength of both fea- 
tures for each of our spectra individually. For each fea- 
ture we fit a Gaussian on a linear continuum background, 
returning the central wavelength. The error is estimated 
by varying the continuum definition either side of the 
feature, performing the fit 100 times, and taking the 
standard deviation of the resulting central wavelength 
distribution. 

The precision is typically ifc2 -5A, adequa t e to d etect 
the shifts of ~10A predicted bv lLentz et al] (|2000f) . We 
found a much greater SN to SN scatter of ±60A, with 
no obvious trend with stretch as might be expected from 
Fig. [ini However, as shown in Fig. [Til for both features 
there is a possible shift with phase in the sense both 
move gradually to longer wavelength as the SN expands. 
Moreover, there is a considerable dispersion at a partic- 
ular phase, particularly before maximum light. 

We find that A2 correlates reasonably well with the 
wavelength of the Si II 4130A - a line commonly used to 
measure the photospheric expansion (Fig. [T^. Thus, to 
first order, the scatter seen in A2 at a particular phase 
appears to arise from variations in kinetic energy within 
the sample. However the phase-dependent trend suggests 
a further process is at work. Fig.[T3]shows the wavelength 
ratio versus phase, coded by stretch, where there is a 
noticeable shift at maximum light (corresponding to a 
redward shift of 50 A in A2). This shift presumably arises 
as a result of the effects of temperature and ionization 
changes which are more pronounced in the UV feature. 

In summary, therefore, the above discussion suggests 
it will be very difficult to use these wavelength features 
as diagnostics of progenitor properties, unless perhaps 
larger samples are available. Clearly the analysis in Lentz 
et al. is oversimplified in several respects. Their models 
considered metallicity variations for a fixed luminosity 
and epoch, whereas it seems a wide range will be needed 
to isolate the effects of metall icity on this part of the 
spectrum. Recent work by iKasen fc WooslevI (|2007[ ) , in 
which a grid of light curve models that span a large range 
in the production of ^^Ni and intermediate mass mate- 
rial (and hence overall kinetic energy), should provide an 
excellent starting place to address these factors in future 
work. 

5.3.3. Host Galaxy and Stretch Dependences 

In view of recent work using the SNLS (|Sullivan et all 
l2006bD where convincing correlations have been found 
between the properties of distant SNe la and their 
host galaxy environment, it is also natural to question 
whether the dispersion implied by the increased scatter 
in Fig.[Tn]is dependent on the stellar population and total 
stellar mass. 

The Keck SN la dataset samples a wide range of passive 
and star-forming host galaxies (Fig. 2]) fo r which stellar 
masses and luminosities are available (see ISullivan et al] 
l2006bl for details). Specific star formation rates and 



masses have been measured for the hosts using SED fits 
to the extr acted colors with the Z-PEG photomet ric red- 
shift code ()Le Borgne fc Rocca-Volmerangell26o^ . 

In Fig. [13] we plot the color-corrected UVl — b color 
against the host-galaxy stellar mass coding the data 
points according to whether the hosts are primarily pas- 
sive or active. The discriminating specific star-formation 
rate was chosen to be 10~^^ Mq yr~^ per unit stellar 
mass (Fig. U). Gas phase metallicities in star forming 
galax ies are we ll-known to correlate with g alaxy mass 
( Tre monti et al.,,2004: Kobulnickv et al.|[2003r ) and recent 
studies have revea led these correlation s were already in 
place at z ~1-1.5 (jShaplev et al.|[2005f) . Over the 3 dex 
range in stellar mass seen in our host galaxies, we would 
expect a change in metallicity of over 1 dex. 

If, as Sullivan et al. surmise, the higher SN la rate in 
active hosts indicate the bulk of those SNe formed rel- 
atively recently, one might expect a tighter or different 
correlation between metallicity and galaxy mass than for 
those in passive galaxies. Unfortunately, Fig. [T3] reveals 
an unavoidable difficulty with the modest Keck sample; 
our four passive galaxies are considerably more massive 
than most of the star-forming galaxies. Thus it is not 
straightforward to separate the effects of galaxy mass and 
specific star formation rate in understanding the implied 
metallicity dispersion. Nonetheless, to the extent we 
can examine any relationships, it seems there are much 
stronger trends between UVl — b and inferred metallic- 
ity for SNe la in star-forming hosts than for those with 
passive hosts. At face value this would imply the most 
metal-rich progenitors have less UV blanketing. Clearly 
a larger sample will be needed to verify these trends. It 
may not be straightforward to sample lower-mass pas- 
sive hosts, given fewer passive hosts are expected at high 
redshift in the hypothesis proposed by Sullivan et al.. 

A more informative comparison of the UV trends can 
be obtained by examining the spectra for those SNe la 
with high and low stretch. Although the host galaxy 
spectral class may contain events wi th a range of stretch 
(c.f. fig. 11 of ISullivan et~al] (|2006bf) V such a categoriza- 
tion of events may be a more direct way of examining 
the spectral properti es of events in old an d young stellar 
populations (see also lBronder et aDl2007t ). 

Fig. [15] shows the mean maximum-light SALT-law 
color-corrected spectra for those events with stretch 
s>l.G3 and s<1.03. We also show the mean spectra split 
by s in the same way when the spectra are color-matched 
(rather than color-corrected) to have the same U — b 
color. These comparisons confirm in detail the photo- 
metric trend noted in the rightmost panel of Fig. [10]- a 
strong and systematic difference in the far UV maximum- 
light spectra of low and high-s SNe la. 

As stretch (or light curve shape) maps closely 
onto both star-forma tion activity in the host galaxy 
([Sullivan et al.l l2006bD and the morphology of the host 
galaxv fe.g. iHamuv et al.ll2000[ ). this comparison pro- 
vides a clear example of the physical differences between 
SNe la in star forming and passive hosts and hence SNe la 
drawn from young and old stellar populations. 

Clearly the mean spectrum of a SN la in a passive 
host with an old stellar population has a depressed far- 
UV continuum compared to one in a more active host 
with a younger stellar population even following color 
corrections based on optical colors. A broadly similar 
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trend is observed by iGuv et al.l (|2007f ). This is in the 
opposite sense than that expected from dust extinction, 
where the higher-stretch SNe la residing in spiral galaxies 
should suffer more extinction and therefore appear red- 
der in the UV. Furthermore, the far-UV colors appear 
more sensitive to stretch than the optical colors. 

Fig. [15] also reveals interesting differences between the 
optical spectral regions of high and low stretch super- 
novae. The majority of these differences originate due 
to the underlying temperature difference, with the high 
stretch supernovae showing stronger high excitation lines 
and ones from doubly-ionized elements. At ~4400 A, a 
line due to Si III 4552 A is found in the high stretch 
supernovae but is absent in the lower stretch examples. 
In addition, the 3700 A Ca-|-Si feature is narrower and 
shows a greater degree of splitting into the separate Ca 
and Si components in the high stretch composite spec- 
trum. Here one can also cl early see the d ifferen ce be- 
tween TZ{Ca II) as defined in lNugent et al.l (|1995| ) where 
the ratio of the emission features on either side of this 
trough are greater in the low stretch composite. Finally, 
the Si II feature at 4130A is shallower i n high-stretch 
SNe c onsistent with trends identified by iBronder et al.l 
(|200l . 

6. DISCUSSION 

One of the original motivations for this study was to 
explore both evolution in the mean properties of SNe la 
over 0< z <0.5 and the dispersion in their UV spectra. 
We have seen remarkably little evidence of evolution in 
the mean spectroscopic properties over this interval, con- 
firming, with higher precision, earlier suggestions based 
largely on lower-S/N data. The primary limitation in 
this aspect of our paper is the paucity of high quality 
UV data at low rcdshift. 

Most of the earlier constraints on the possible evolu- 
tion of SNe la were largely based o n considera t ion o f 
photometric measures. For example, iRiess et aD ()2004l ) 
derived the rest-frame U — B and B — V distributions of 
a heterogeneous sample of SNe la at maximum light over 
z<1.5 as a proxy for the more precise constraints possi- 
ble with spectroscopic data. They ruled out any color 
evolution greater than 0.02 mag in the mean in U — B. 

Spectroscopic data is substantially more precise in its 
tracking of evolution and dispersion than broad-band 
photometry for two reasons. Foremost, detailed dif- 
ferences in key line diagnostics are lost in photometric 
data. Secondly, no assumptions need be made about k- 
corrections in gathering and comparing data over a broad 
redshift range. 

At the time this Keck survey was underway, there was 
no comprehensive spectroscopic dataset against which 
our UV data could b e compared, other t han the mean 
template published bv lNugent et al.l ()2002D and discussed 
in § 14.31 However, very recently, as a result of a survey of 
high-redshift SNe la with the ACS g rism onboard Hubble 
Space Telescope, iRiess et all (|2007D have published the 
mean rest-frame UV spectrum of 13 SNe la more distant 
(z > 1) than those studied here; the average redshift 
of this sample is z=1.3. In this study, used to demon- 
strate the presence of dark energy at 2; >1, it is claimed 
that the sample-averaged spectral energy distributions 
observed at z >1 are consistent with that observed lo- 
cally and that any spectral evolution is still undetected. 



However, no quantitative statement is made to support 
this claim. 

We believe t h at th e spectral comparison undertaken 
by IRiess et al.l ()2007D is less precise than that under- 
taken in this study, and hence less valuable as a means 
of justifying the continued use of local relations in con- 
structing SN la Hubble diagrams. Foremost, the ACS 
grism data samples SNe la over a much wider range of 
phase than the comparisons undertaken here. Although 
the individual spectra comprisi ng the mean ^=1. 3 ACS 
spectrum are not tabulated by .Riess et al.l (|2007t ) , their 
phases range over at least 15-20 days and almost all are 
post-maximum light, compared with the narrow phase 
range maximum-light (-4 to -1-4 days) and pre-maximum 
light (<-4 days) comparisons presented here. A second 
consideration is the inevitable poor signal-to-noise ratio 
of these very high redshift spectra which precludes de- 
tailed comparisons. 

As we have shown in F igs. |5| and |6l the local tem- 
plate (jNugent et al.ll2002t l is a poor basis from which to 
make such comparisons. More representative datasets 
are needed for reliable claims. Internally within our own 
data. Fig. |7| reveals little evolution although the redshift 
baseline is small. It will be important to examine the 
case for evolution from 0.2< z <1.5 by combining in a 
consistent manner both the IRiess et al.l (2007) dataset 
with that presented here, as well as with higher-redshift 
SNLS SNe la (out to z ~ 1) observ ed during the routin e 
survey spectroscopic screening (e.g. IBronder et al.ll2007[ ) . 

The second significant finding in this study is the 
marked increase in the dispersion among our SNe la 
shortward of 3300A (Fig. [T0|) . Although the scatter in 
U is comparable to that seen in local data, it increases 
significantly at shorter wavelengths. Even allowing for 
differing amounts of dust extinction within each host 
galaxy, this amounts to more than a factor of 2 vari- 
ation in continuum flux at maximum light. Although 
theoretical models predict a strong sensitivity to metal- 
licity variations at this wavelength, our variations are 
also considerably larger than those predicted. We have 
demonstrated that a color correction based on the B ~V 
color of the SN and either a CCM Milky Way extinction 
law or a SALT color law, can only marginally reduce, 
and not eliminate, this UV dispersion. 

We have found it hard to isolate the physical causes 
of this significant UV dispersion. We confirm earlier 
work at longer wavelengths (A ~3500-4000A) that shows 
that stretch (or equivalently host galaxy class) is par- 
tially responsible. However at the shortest wavelengths 
(A ~3000-3300A corresponding to our UVl diagnostic), 
additional effects are clearly important that are not ac- 
counted for by the color-correction techniques in use in 
current cosmological programs. In addition to the newly- 
found intrinsic scatter at short wavelengths, new trends 
with phase are also seen in the wavelengths of diagnostic 
features in this region. 

Redward of 4000 A , corrections for stretch or color 
work well in normalizing SNe la, however blueward of 
4000A, significant scatter remains even after such correc- 
tions are made. This may be related to a change in the 
dominant source of opacity in SNe la. Redward of 3500- 
4000 A electron scattering opacity dominates, but at UV 
wavelengths a forest of overlapp ing lines is the dominant 
source of opacity (see Fig. 1 of iHillebrandt fc Niemeveil 
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(|2000f )). Electron scattering is a continuous process in- 
volving well-understood physics, but line opacity de- 
pends sensitively on abundances, ionization states, and 
possibly non-LTE effects. 

Fig. [16] illustrates this point by comparing the 
wavelength dependence of the line and electron scat- 
tering opacity a t max imum light for a model by 
iKasen k. WooslevI ()2007t ) that provides a good match for 
a normal SN la (Kason, priv. commun.). The data in 
question refers to that at a depth of 7000 kms~^ where 
intermediate mass and Fe peak material are well-mixed. 
A common feature in these models is the drop in line 
opacity compared to the electron scattering opacity near 
4000 A as seen here. Thus it is understandable that 
the emerging UV flux is highly susceptible to changes in 
the line opacity (due to initial conditions and/or mate- 
rial synthesized during the explosion) while the optical 
and near-IR spectral behavior are dominated by electron 
scattering opacity at this phase. 

What are the possible consequences of the above vari- 
ations in terms of the use of SNe la as probes of the ex- 
pansion history? In the highest redshift surveys, includ- 
ing those proposed with future facilities, cross-color k- 
corrections are needed to estimate rest-frame light curves 
from the observations, typically undertaken in the far-red 
and near-infrared for z >1.5 (Ricss ct al. 2004, 200'^. 
This necessitates the adoption of a suitable template 
whose SED is reliable in the 3000-4000A region. At the 
most fundamental level, the dispersion in our UV spec- 
tra will contribute a statistical uncertainty at the 0.05- 
0.1 magnitude level depending on how precisely particu- 
lar observed filters match to the chosen rest-frame band- 
passes. 

A more worrying bias would follow the adoption of an 
incorrect template. Possibly the most significant finding 
in our work is the demonstration that the mean UV spec- 
trum is different for high and low stretch SNe la (Fig. [T5)) 
and thus, presumably, for those that occur in passiv e 
and actively star- forming galaxies (jSullivan et al.|[20"06bD . 
Adopting a single template would lead to a systematic 
bias which could become increasingly serious at high red- 
shift where the mix between the two populations changes 
and high-stretch SNe, typically originating from younger 
progen itor systems, becom e increasingly common. Re- 
cently, iHowell et al.l (|2007( ) have demonstrated the ex- 
istence of an increased fraction of high-stretch SNe la 
at higher redshifts above and beyond the selection effect 
that high stretch, brighter SNe la are easier to detect. 

Using our individual maximum light spectra, we can 
evaluate the impact of the UV dispersion on sample 
cross-color /c-corrections, addressing both the statistical 
error arising from the intrinsic scatter and the systematic 
difference arising from the stretch-based spectral differ- 
ences we have found. Fig. [T71 shows the result. 

For a typical future dark energy experiment, based on 
securing the equation of state parameter w to 5% us- 
ing z >1 SNe la, photometric corrections of better than 
±0.02 mag are required This requirement is indicated 
on Fig. [T71 

In terms of the statistical error, the observed dispersion 
is typically ±0.05-0.10 mag, larger than that required by 

e.g. White Papers submitted to the US Dark Ene rgy Task 
Force Committee; |http : / / www. nsf . gov/mps /ast/detf.jsp| 



a factor of several. However, such a dispersion, if ran- 
domly distributed, need not present a fundamental ob- 
stacle to progress with a survey spanning a wide redshift 
range utilizing a large number of SNe la. As Fig. [17] 
shows, at certain redshifts where the rest-frame and ob- 
served filters closely match, the dispersion has negligible 
effect and SNe la at these redshifts could be more heavily 
weighted. 

The more worrying trend, parti cularly given t he rec ent 
demonstration of stretch bias bv lHowell et al.l ()2007[ ). is 
the likelihood of a systematic error introduced by adopt- 
ing an incorrect template for the redshifts where the rest- 
frame and observed filters do not overlap. Fig. [TTl shows 
the potential of this error via a comparison of the disper- 
sion independently for high and low stretch events. The 
differences between these two categories are comparably 
large and indicate the importance of securing a physical 
understanding of the UV variations seen in our survey. 

7. CONCLUSIONS 
We summarize our findings as follows: 

1. We have secured high signal-to-noise ratio Keck 
spectra for a sample of 36 intermediate redshift 
SNe la, observed at various phases, spanning the 
redshift range 0.15< z <0.7, and drawn from the 
Supernova Legacy Survey (SNLS). We demonstrate 
via inspection of the SN properties that our Keck 
sample is a reasonably fair subset of the larger sam- 
ple of distant SNe la being studied by the SNLS. 

2. We develop a new method for removing host galaxy 
contamination from our spectra based on measures 
of the galaxy and SN photometry. These refine- 
ments to traditional spectral reduction techniques 
allow us to achieve host-galaxy subtracted and flux- 
calibrated rest-frame spectra of high quality, ex- 
tending down to rest-frame wavelengths of 2900A. 

3. Although no strong evidence is found for spectral 
evolution in the mean early-phase and maximum 
light spectra, when compared to local data, such 
evolutionary tests are hampered by the paucity of 
quality data at low redshift and a significant scat- 
ter in the spectra shortward of 4000A. We find no 
evidence for evolution internal to our data. We ar- 
gue that the well-us ed local UV spectral template 
([Nugent et al.l[200^ is likely to be less represen- 
tative than the mean spectrum compiled from the 
Keck data which we tabulate with the measured 
dispersion for use in future cosmological applica- 
tions. 

4. Our principal finding is a large scatter from one 
SN to the next in the rest-frame UV spectrum 
even after standard dust corrections are made. 
By constructing various photometric bandpasses 
that avoid uncertainties arising from differential k- 
corrections associated with the range of redshifts 
in our sample, we show that while we can repro- 
duce the stretch-dependent trends seen locally at 
3500-4000A, the scatter at 3000-3400A is 3-5 times 
larger. 
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5. Although progenitor metallicity may drive some of 
the trends seen in the Keck data, the UV varia- 
tions are much larger than in contemporary models 
which span the expected metallicity range. More- 
over, the UV spectrum also changes with phase in 
a manner which is not consistent with models. We 
conclude there are significant variations in the UV 
properties of SNe la which are not accounted for 
by either the presently-employed empirical trends 
or the available SN la models. 



6. As an illustration of the importance of understand- 
ing these new results, we calculate the error aris- 
ing from the use of a single UV spectral template 
for calculating the cross-color k correction, a cor- 
rection essential for constructing the SN la Hub- 
ble diagram as a probe of the expansion history. 
The dispersion arising from our UV spectra, if not 
randomly distributed along the Hubble diagram, 
presents an uncertainty 2-3 times larger than than 
would be necessary for recovering the equation of 
state parameter w to 5% using SNe la at z ~1. We 
conclude that further detailed studies are essential 
if SNe la are to be useful for precision measures of 
dark energy. 
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TABLE 1 

Distant Supernova Sample 



SN Namp" 


RA 


DEC 






O Li L LLli 


iNcaicoii y 


/o 








t Tonnn\ 
(JzUOUJ 




(Days) 




pnot (UayH ) 


Inc. 


Source 


OSDlau 


02:24:10.38 


-04:02:14.9 


0.5043 ± 0.0005 


-1.6 + 0.3 


1.14 + 0.02 


+0.6 


432 


Galaxy 


OSDlaw 


02:24:14.78 


-04:31:01.6 


0.582 ± 0.005 


+2.4 + 0.3 


1.08 + 0.02 


+0.6 


856 


Galaxy 


OSDlco 


02:26:16.23 


-04:56:05.7 


0.679 + 0.001 


+7.1 + 0.5 


1.07 + 0.04 


+0.7 


321 


Galaxy 


OSDldj 


02:26:19.08 


-04:07:09.3 


0.40 + 0.01 


-1.8 + 2.3 


1.25 + 0.61 


+0.1 


4199 


Template Fit 


03D3af 


14:21:14.92 


+52:32:15.3 


0.5320 + 0.0005 


+2.8 + 0.4 


1.02 + 0.03 


+3.8 


368 


Galaxy 


03D3aw<= 


14:20:53.61 


+52:36:20.6 


0.449 + 0.001 


-0.8 + 0.2 


0.94 + 0.03 


-13.7 


52 


Galaxy 


03D3ay<= 


14:17:58.43 


+52:28:57.4 


0.3709 + 0.0003 


-1.0 + 0.2 


1.00 + 0.03 


-14.6 


72 


Galaxy 


03D3ba<= 


14:16:33.44 


+52:20:32.1 


0.2912 + 0.0003 


+12.8 + 0.4 


1.04 + 0.03 


-15.5 


259 


Galaxy 


03D3bbf 


14:16:18.78 


+52:14:55.3 


0.2437 + 0.0003 


+2.2 + 0.3 


1.18 + 0.03 


-16.2 


2323 


Galaxy 


03D3bh<= 


14:21:35.89 


+52:31:37.6 


0.2486 + 0.0003 


-3.5 + 0.2 


1.02 + 0.03 


-16.1 


140 


Galaxy 


03D3bl 


14:19:55.90 


+53:05:51.0 


0.3553 + 0.0005 


+3.7 + 0.5 


0.98 + 0.03 


-0.1 


270 


Galaxy 


03D3CC 


14:19:45.25 


+52:32:25.3 


0.4627 + 0.0003 


+8.5 + 0.4 


1.04 + 0.02 


-0.1 


204 


Galaxy 


03D3cd 


14:18:39.95 


+52:36:43.8 


0.4607 + 0.0005 


-5.6 + 0.1 


1.15 + 0.02 


-0.1 


625 


Galaxy 


03D4ag 


22:14:45.79 


-17:44:23.0 


0.2847 + 0.0003 


-5.0 + 0.1 


1.08 + 0.02 


-1.6 


141 


Galaxy 


03D4cjf 


22:16:06.66 


-17:42:16.7 


0.27 + 0.01 


-7.4 + 0.1 


1.09 + 0.01 


+1.5 


5000 


Template Fit 


03D4dh 


22:17:31.04 


-17:37:46.9 


0.6268 + 0.0004 


+0.5 + 0.3 


1.08 + 0.02 


+0.6 


444 


Galaxy 


03D4gl'= 


22:14:44.17 


-17:31:44.4 


0.571 + 0.001 


+7.7 + 0.5 


1.00 + 0.03 


+14.6 


347 


Galaxy 


04Dlhd 


02:26:08.85 


-04:06:35.2 


0.3688 + 0.0004 


+1.6 + 0.1 


1.06 + 0.01 


+1.4 


1998 


Galaxy 


04Dljg 


02:26:12.56 


-04:08:05.3 


0.5842 + 0.0005 


-2.0 + 0.2 


1.05 + 0.02 


+1.3 


128 


Galaxy 


04Dloh<i 


02:25:02.37 


-04:14:10.5 


0.59 + 0.01 


-6.0 + 0.2 


1.03 + 0.02 


+0.1 


512 


Template Fit 


04Dlrh 


02:27:47.16 


-04:15:13.6 


0.4349 + 0.0005 


+2.6 + 0.3 


1.02 + 0.03 


-0.1 


393 


Galaxy 


04Dlsk 


02:24:22.77 


-04:21:13.3 


0.6634 + 0.0005 


-0.4 + 0.6 


0.86 + 0.04 


+0.0 


139 


Galaxy 


04D2gc 


10:01:39.26 


+01:52:59.5 


0.5216 + 0.0005 


-0.2 + 0.3 


1.14 + 0.03 


-0.9 


367 


Galaxy 


04D2kr 


10:00:37.32 


+01:42:43.1 


0.7441 + 0.0005 


-2.5 + 0.6 


1.06 + 0.04 


+0.6 


91 


Galaxy 


04D3cp<^ 


14:20:23.95 


+52:49:15.5 


0.83 + 0.02 


+5.2 + 0.4 


1.07 + 0.02 


+0.5 


557 


Galaxy 


04D3ez 


14:19:07.91 


+53:04:18.8 


0.2630 + 0.0005 


+1.6 + 0.1 


0.92 + 0.02 


-0.4 


78 


Galaxy 


04D3fk 


14:18:26.21 


+52:31:42.7 


0.3578 + 0.0007 


-6.9 + 0.1 


0.97 + 0.01 


-0.4 


136 


Galaxy 


04D4in 


22:15:08.58 


-17:15:39.8 


0.5160 + 0.0005 


-4.6 + 0.1 


1.13 + 0.01 


+2.0 


443 


Galaxy 


04D4jr 


22:14:14.33 


-17:21:00.9 


0.482 + 0.007 


-0.2 + 0.2 


1.16 + 0.02 


+2.7 


5000 


Template Fit 


05Dlhk' 


02:24:39.16 


-04:38:03.0 


0.2631 + 0.0002 


-8.6 + 0.1 


1.16 + 0.01 


-0.8 


337 


Galaxy 


05Dlhn 


02:24:36.25 


-04:10:54.9 


0.1489 + 0.0006 


-7.2 + 0.1 


1.06 + 0.02 


-0.8 


232 


Galaxy 


05Dlif 


02:24:29.71 


-04:34:13.0 


0.763 + 0.001 


-5.0 + 0.3 


1.03 + 0.03 


-0.0 


5000 


Galaxy 


05Dlix 


02:24:19.95 


-04:40:11.7 


0.49 + 0.01 


-8.7 + 0.2 


1.05 + 0.01 


-0.7 


5000 


Template Fit 


05Dliy 


02:27:39.97 


-04:25:21.3 


0.2478 + 0.0003 


-9.2 + 0.1 


0.88 + 0.01 


-0.7 


52 


Galaxy 


05D21e 


10:01:54.85 


+02:05:34.7 


0.7002 + 0.0005 


+6.3 + 0.4 


1.10 + 0.04 


-0.0 


2647 


Galaxy 


05D2mp 


09:59:08.61 


+02:12:14.6 


0.3537 + 0.0004 


-3.9 + 0.2 


1.11 + 0.02 


+0.7 


831 


Galaxy 



^ The SNLS ID format is explained in iHowell et all ( |2005l) ; briefly, it relates to the year of discovery (03, 04 etc.), the relevant 
SNLS field (Dl through D4), and a running ID (e.g. aa, ab, ac ctc.).^ "Effective" rest-frame phase relative to maximum light i.e. 
(T'oba — 7\nax) / {s x z))^ Estimated through the spectral extraction aperture; increases greater than 5000% are listed as 5000%.^ Light 
curve deemed inadequate in temporal coverage. ° Spectrum deemed inadequate in signal-to-noise ratio. ^ Spectroscopically peculiar; 
excluded from primary analysis sample. 
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TABLE 2 
Observing log 



SN Name Date Obs 


MJD 


Seeing 


Slit 


Dichroic Exptime SN apparent 






(") 


(") 


(s) i' mag. 



03Dlau 


2UU3-Uy-23 


52905 


.51 


0. 


03Dlaw 


2UU3-Uy-23 


52905 


.58 


0. 


03Dlco 


ZOuo- 11-22 


52965 


.48 


0, 


03Dldj 




52964, 


.47 


0, 


03D3af 


2UU3-U4-U7 


52736 


.41 


0, 


03D3aw 


2UU3-U5-UD 


52765 


.53 


0, 


03D3ay 


2UU3-U5-UD 


52765 


.47 


0, 


03D3ba 


20U3-U5-UD 


52765 


.38 


0, 


n o o 1, 

UoDJbb 


2UU3-U5-UD 


52765 


.29 


0, 


03D3bli 


2()U3-U5-Uo 


52765 


.32 


0, 


03D3bl 


2UU3-UD-U1 


52791 


.33 


0, 


03D3cc 


2UU3-UD-U2 


52792, 


.27 


1. 


03D3cd 


2003-06-01 


52791, 


.27 


1. 


03D4ag 


2003-07-02 


52822, 


.54 


0. 


03D4cj 


2003-08-27 


52878, 


.32 


0. 


03D4an 


2003-09-23 


52905, 


.43 


0. 


03D4gl 


nr\no 1 1 oo 

2003-11-22 


52965 


.26 


0. 


04Dlhd 


2004-09-21 


53269 


.47 


0, 


04Dljg 


2004-09-21 


53269 


.56 


0, 


U4Dlon 


onn /I 1 n 1 ri 

zUU4-iU-i9 


53297 


.54 


1, 


04Dlrh 


i^r\f\ A 1 O 1 ■! 

2004-12-14 


53353 


.24 


0, 


UrtJ— ' J-oiS. 


9004 19 14 


ooooo 


.36 


0, 


04D2gc 


2004-04-22 


53117, 


.00 


0, 


04D2kr 


2004-12-14 


53353, 


.58 


0, 


04D3CI) 


2004-04-23 


53118, 


.45 


0, 


04D3ez 


2004-04-22 


53117, 


.00 


0, 


04D3fk 


2004-04-22 


53117, 


.00 


0, 


04D4in 


2004-09-21 


53269, 


.31 


0, 


04D4jr 


2004^10-19 


53297, 


.28 


1. 


OSDlhk 


2005-11-30 


53704, 


.37 


0. 


OSDlhn 


2005-11-30 


53704, 


.40 


0. 


05Dlif 


2005-12-01 


53705, 


.31 


0. 


OSDlix 


2005-11-30 


53704, 


.28 


0. 


05Dliy 


2005-11-30 


53704, 


.43 


0. 


05D21e 


2005-12-01 


53705, 


.57 


0. 


05D2mp 


2005-11-30 


53704, 


.53 


0. 



8 


1.0 


560 


5400 


22.4 


8 


1.0 


560 


5400 


22.7 


9 


1.0 


680 


9000 


23.5 


8 


1.0 


560 


6000 


22.0 


9 


1.0 


680 


5400 


22.7 


9 


1.0 


680 


4800 


22.2 


9 


1.0 


680 


3600 


21.7 


9 


1.0 


680 


3600 


21.7 


9 


1.0 


560 


2000 


19.8 


9 


1.0 


560 


3600 


20.9 


9 


1.0 


680 


3600 


22.0 





1.0 


680 


5400 


22.1 





1.0 


680 


9777 


22.2 


9 


1.0 


560 


5400 


21.1 


9 


1.0 


560 


4000 


21.2 


8 


1.0 


560 


3600 


22.6 


9 


1.0 


680 


7200 


22.6 


9 


1.0 


560 


6000 


21.6 


9 


1.0 


560 


6000 


22.5 





1.0 


560 


7200 


22.9 


9 


0.7 


560 


6000 


21.9 


9 


1.0 


560 


9000 


23.2 


7 


0.7 


560 


6600 


22.5 


9 


1.0 


560 


9600 


22.8 


8 


0.7 


560 


2700 


23.1 


8 


0.7 


560 


3600 


21.0 


8 


0.7 


560 


4500 


22.3 


9 


1.0 


560 


7200 


22.4 





1.0 


560 


7200 


22.0 


9 


1.0 


560 


2400 


21.5 


9 


1.0 


560 


2100 


20.8 


9 


1.0 


560 


9000 


23.2 


9 


1.0 


560 


7200 


22.9 


9 


1.0 


560 


4800 


22.1 


9 


1.0 


560 


7200 


23.2 


9 


1.0 


560 


6600 


21.9 
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TABLE 3 

Mean SN Ia spectrum (Early) 



Wavelength (A) 


Flux 


Flux lower 


Flux upper 


2805 


0.506 


0.302 


0.629 


2810 


0.579 


0.322 


0.643 


2815 


0.417 


0.364 


0.662 


2820 


0.477 


0.390 


0.682 


2825 


0.482 


0.431 


0.708 


2830 


0.529 


0.449 


0.724 


2835 


0.747 


0.472 


0.757 


2840 


0.723 


0.483 


0.792 


2845 


0.676 


0.481 


0.813 


2850 


0.719 


0.495 


0.830 


Note. — The complete version 


of this table 


is in the electronic edition of the Journal. The 


printed edition contains only a sample. 




TABLE 4 




Mean SN Ia spectrum (Max) 


Wavelength (A) 


Flux 


Flux lower 


Flux upper 


2805 


0.287 


0.262 


0.411 


2810 


0.360 


0.282 


0.422 


2815 


0.372 


0.299 


0.435 


2820 


0.438 


0.317 


0.448 


2825 


0.355 


0.334 


0.456 


2830 


0.477 


0.349 


0.470 


2835 


0.433 


0.357 


0.472 


2840 


0.378 


0.367 


0.483 


2845 


0.442 


0.375 


0.496 


2850 


0.428 


0.377 


0.506 



Note. — The complete version of this table 
is in the electronic edition of the Journal. The 
printed edition contains only a sample. 
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Fig. 1. — Redshift (left) and phase (right) distributions of the high-redshift SNLS SNe la studied in this paper. The vertical dashed 
lines in the phase distributions show how the SN la sample is divided into "early" and "maximum-light" spectra. Numbers (and the open 
histograms) refer to the primary sample of 36 high signal-to-noise ratio spectra whereas those in parentheses (and the shaded histograms) 
refer to the sample of 26 used to construct the mean UV spectra (see text). 



5500 6000 



Original data 

1 .50 — HosI galaxy 

Final color-adjusted speotr 




Rest Wavelength (A) 
3000 3500 4000 4500 



5000 5500 6000 6500 
Observed Wavelength (A) 




5500 6000 6500. 
Observed Wavelength (A) 



Fig. 2. — Examples of the host galaxy subtraction and flux calibration techniques for two SNe la suffering host contamination. The blue 
spectrum represents the observed data, the red spectrum the estimated host galaxy spectrum from fits to the broad-band galaxy u*g'r'i'z' 
photometry (red circles), and the black spectrum is the host-galaxy subtracted final spectrum color-corrected to the g'r'i'z' SN photometry 
(black circles). A smoothed version of the final spectrum is over-plotted, the heavy line denoting the spectral region deemed to have the 
most reliable flux calibration. The photometric points are plotted at the effective wavelengths of the fllters on the appropriate spectrum, 
and the total system response through the g'r'i'z' fllters overplotted at the observed wavelengths. The position of the dichroic between 
the two arms of the spectrograph is also marked. Left (a): SNLS-04D3ez a,t z = 0.263. This SN had an increase in i' of ~ 80% on the day 
of observation. Right (b): SNLS-03D3cc at z = 0.463, with an increase in i' of ~ 200%. Note that due to the color correction applied to 
the final spectrum, as plotted the black spectrum differs slightly from the simple subtraction of the red from the blue. [See the electronic 
edition of the journal for Figs. [2]3-|2]26] 



Verifying the Cosmological Utility of SNe la 



19 
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Spec. Norm, (b) 
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Fig. 3. — Maximum light SN la spectra from lLentz et all (j2000') at a variety of different metallicities (top to bottom: 0.1, 0.3, 1.0, 3.0 and 
10.0 times solar metallicity). Over-plotted are the standard Bcssell (1990) U and B filter responses, as well as the UV box filters ("UVl" 
and "UV2") and spectroscopic normalizing filter ("b") used in the analysis. Features thought to represent possible metallicity diagnostics 
(Ai, A2) are also marked, as well as the position of the blue-shifted Si II 4130A absorption feature. 




Fig. 4. — Comparison of SN stretch, SN B-V color, host galaxy specific star formation rate and SN/galaxy projected separation for the 
Keck SN la sample (black histogram; present paper) and the entire SNLS sample within the appropriate redshift range discovered over the 
same period (grey histogram). Galaxies with undetected star formation were placed at 10"-*^^ Mq yr"-"^ per unit stellar mass. 
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2.0 




0.0 1 , , I , , , , I , , , , I , , , , I , , , , I , , , , I 

3000 3500 4000 4500 5000 5500 

Wavelength (A) 

Fig. 5. — The mean hig h-redshift maximum light (effective day < ± 4 days) rest-frame UV SN la spectrum compared to the local 
average template of lNugent et al. ( 2002). Over-plotted in light grey are 100 bootstrap-resampled mean spectra drawn from the high-redshift 
population; the dotted lines show the region containing 90% of this distribution. The local template has been color-adjusted to match the 
high redshift data. 
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Fig. 6. — As Fig. [5] but showing the pre-maximum (effective day <-4 days) high-redshift spectrum versus the local template. 
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Fig. 7. — The mean high-redshift maximum light rest-frame UV SN spectrum split into two rcdshift intervals below (black) and above 
(red) z=0.5. 
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Fig. 8. — A test of the validity of applying a color correction to those Keck spectra sampled at maximum light. The top panel shows the 
mean observed spectrum for two subsamples split according to the rest-frame B — V color at maximum light; red corresponds to SNc with 
B — V > —0.057, blue to those with B — V <-0.057. The middle panel shows the same comparison after applying a Cardelli ct al. (1983) 
reddening correction to each individual spectrum using an E{B — V) estimated from the light curve fits, and Rg = 4.1. Varying Rg has 
a negligible effect on this comparison. The lower panel shows the comparison after correcting individual spectra using the SALT color law 
IIGuv et al.ll2005l) . 
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Fig. 9. — The deviation of the 15 maximum- Ught spectra from the mean spectrum in cr-units derived from their error spectra 15.2.21 
The top panel shows the uncorrected spectra, the middle panel the spectra corrected with the CCM dust law, and the lower panel spectra 
corrected with the SALT color law. The horizontal lines denote ±2-(t deviations. 
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Fig. 10. — Various SN la stretch-color relations, with colors measured directly from the SALT color law corrected Keck spectra. The 
colors shown are U — b, UVl — b and UV2 — fe - see Fig.[3]for filter definitions. SNe are color-coded according to the phase of the spectrum 
as defined in Fig . [T] The black p oints in the U — b distribution refers to the U — b distribution of local SNe la at maximum light from local 
SN surveys fe.e. lJha et al.ll2006l) . The error-bars for the SNLS SNe include propagated uncertainties from host galaxy subtraction 13.31 1. 
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Fig. 11. — Variation of wavelength with phase for the two UV diag nostic features mar ked in Fig.|3] a clear trend is apparent. The vertical 
line refers to the full extent of the metallicity variation predicted by ILentz et al.l 112000) . 
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Fig. 12. — Comparison between the wavelength of the Si II 4130A feature commonly used as a measure of the photospheric expansion 
velocity and that of the UV diagnostic A2 (see Fig. [3] for definition). Spectra taken before maximum light are indicated with blue symbols; 
those after maximum light with red symbols. Numbers indicate the phase in rest-frame days. 
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Fig. 13. — Ratio of the wavelength of the Si II 4130A feature to the UV diagnostic A2 versus phase with data points sized according to 
the stretch. There is a clear discontinuity in the behavior after maximum light. 
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Fig. 14. — The variation of the Keck SN la UVl — b spectral color as a function of the stellar mass of the SN host galaxy color coded 
according to the specific star formation rate. The division is made at a rate of 10~^^Mq yr~^ per unit stellar mass. The color error-bars 
include propagated uncertainties from host galaxy subtraction. The conversion from stellar mass to gas-phase metallicity using the relation 
of .Trcmonti ct al.. (.2004 ) is also shown. 
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Fig. 15. — The mean SALT-law color-corrected (left) and color-matched (right) Keck spectra for SNe la split according to the SN light 
curve stretch, s. The red spectrum refers to events with s < 1.03, the blue spectrum to those with s > 1.03. 
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Fig. 16. — Comparison of the wavelength-dependent line and electron scattering opacity for a typical model in lKaseii fc WooslevI 120071) 
at peak SNIa brightness. The data refers to a depth of 7000 kms~^. Note the drop in line opacity with respect to the electron scattering 
opacity near 4000 A . This behavior makes the emergent UV flux highly sensitive to changes in the line opacity whereas the optical and 
near-IR spectral regions are largely dominated by electron scattering opacity. 
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Fig. 17. — Redshift dependence of various maximum light cross-filter fc-corrections differenced to that based on the local template. Lines 
represent deviations observed for individual Keck spectra color-coded according to stretch: red lines refer to events with s < 1.03, while 
blue to those with s > 1.03. The dashed lines indicates the approximate precision necessary to secure an equation of state parameter w to 
a precision of 5% using z >1 SNela (see text for discussion). 



